Ab Initio Theory of Scattering-Independent Anomalous Hall Effect 
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We report on first-principles calculations of the side-jump contribution to the anomalous Hall 
conductivity (AHC) directly from the electronic structure of a perfect crystal. We implemented 
our approach for a short-range scattering disorder model within the density functional theory and 
computed the full scattering-independent AHC in elemental bcc Fe, hep Co, fee Ni, and Llo FePd 
and FePt alloys. The full AHC thus calculated agrees systematically with experiment to a degree 
unattainable so far, correctly capturing the previously missing elements of side-jump contributions, 
hence paving the way to a truly predictive theory of the anomalous Hall effect and turning it from 
a characterization tool to a probing tool of multi-band complex electronic band structures. 

PACS numbers: 72.25.Ba, 72.15.Eb, 71.70.Ej 



The anomalous Hall effect (AHE) in ferromagnets is 
one of the most celebrated transport phenomena in solid- 
state physics [1]. It has been researched intensely in the 
past decade after it was realized that the intrinsic con- 
tribution (IC) could be interpreted in terms of the Berry 
phases of Bloch electrons in a solid [H, For almost 
ten years the IC was the only one, which could be ac- 
cessed in density functional theory (DFT) calculations 
of the AHE The ability to estimate the impurity- 

driven (i.e. extrinsic) contributions to the AHE has been 
remarkably limited so far, thus hindering the predictive 
power in understanding and engineering the AHE trans- 
port properties of real materials. 

Besides the IC, a^"*, there are two extrinsic disorder- 
driven contributions to the anomalous Hall conductivity 
(AHC). In metallic systems they can be distinguished ac- 
cording to their parametric dependency on the impurity 
concentration n, with skew-scattering conductivity, a®^, 
proportional to 1/n Q, and the side- jump contribution 
(SJC), cr^j, independent of impurity concentration [H. 
The fact that the SJC, although originating from im- 
purities, does not depend on their concentration, makes 
it one of the most challenging electron scattering mecha- 
nisms to understand and suggests a close relation to the 
intrinsic contribution of the AHE. This behavior arises in 
metallic systems from the I/cft expansion of the trans- 
port coefficients in Hnear response for coupled multi-band 
systems. Consequently, since the SJC does not depend 
on the disorder strength, it cannot be easily separated 
from the IC in low temperature dc measurements Q. 
So far, only for an Llo-ordered FePd ferromagnetic al- 
loy clear evidence has been presented that the side-jump 
contribution can dominate over other mechanisms of the 
AHE in a wide temperature range through a comparison 
of theoretical IC and experimentally measured AHC 

Previous numerical ab initio treatments have mainly 
concentrated on the IC Q and defined the SJC by tak- 
ing the zero disorder limit of coherent potential approxi- 



mation disordered alloys calculations p^] . Within such a 
treatment the indirect computation of the SJC presents a 
significant computational challenge considering also that 
the exact knowledge of the disorder potential in the sys- 
tem is necessary in this case. Since very often the exper- 
imental data are obtained on samples with unknown im- 
purity content and disorder type, it is highly desirable to 
be able to evaluate the SJC expHcitly from the electronic 
structure of a perfect crystal. A direct computation of 
the SJC in models with disorder was not feasible until 
recently when it was shown that, assuming short-range 
uncorrelated disorder model, the SJC may be indeed cal- 
culated directly from the ideal electronic structure of a 
crystal without any disorder [11|, rendering possible a 
rigorous numerical study that can be fully compared to 
experimental results and can set the stage to a truly pre- 
dictive theory of the AHE. While the validity of derived 
expressions has been demonstrated for simple models, 
values for the scattering-independent SJC in fundamen- 
tal ferromagnetic materials such as Fe, Co or Ni, have 
not been obtained so far. 

In this Letter we report on calculations from first prin- 
ciples of the values of the scattering-independent side- 
jump conductivity in elemental bcc Fe, hep Co, fee Ni, 
as well as ordered FePd and FePt alloys directly from 
the electronic structure of their pristine crystals. Our 
calculations unambiguously show that the calculated val- 
ues of the total scattering-independent contributions, IC 
and SJC, agree systematically with experiments to a level 
that was not reached before in these materials. More 
importantly, the inclusion of the scattering-independent 
SJC accounts consistently for the discrepancy between 
the IC and the measured values in a very non-trivial fash- 
ion. We analyze the side-jump as a Fermi surface prop- 
erty and demonstrate that it shows a strong anisotropy 
with respect to the magnetization direction in the crystal, 
even more pronounced than that found for the intrinsic 
AHC 



12|, Il3| . Additionally, we demonstrate the impor- 
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tance of correlation effects in describing the AHE in fee 
Ni correctly. 

The starting point for the theory of the scattering- 
independent side-jump are the retarded Green's func- 
tion in equilibrium and the Hamiltonian of a general 
multiband noninteracting system in three spatial dimen- 
sions. At the first step we expand the self-energy of the 
system Seq in powers of potential V^(r), which describes 
scattering at impurities. For a short-range scattering dis- 
order model, scalar delta-correlated Gaussian disorder or 
delta-scattering uncorrelated disorder, the contribution 
to the self-energy which is of first order in V^(r) vanishes, 
since one can assume that {V{r)) = or else absorb 
(V(r)) into the Hamiltonian, a procedure which results 
in a simple shift of the energy levels. Further, inserting 
the expression for the self-energy within these simple dis- 
order models into appropriate equations for the current 
densities derived following the Kubo-Streda formalism, 
rotating into eigenstate representation and keeping only 
the leading order terms in the limit of vanishing disorder 
parameter V, i.e. ignoring skew-scattering contributions, 
the scattering-independent part of the AHE conductivity 
may be written as a^^^ = a^"* + a^J, where 



^lm}2{fn-fm)- 



can be recovered as the intrinsic contribution In 
this expression indices n and m run over all bands with 
occupations fn and respectively, Vnm,i are the ma- 
trix elements of the velocity operator Vi = dhkiH^ and 
uoniy) = £n(k)/^ with SnO^) as band energies. The 
scattering-independent SJC to conductivity cr*^^^ reads for 
inversion-symmetric systems: 
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Here £f is the Fermi energy and the imaginary part of the 
self-energy ImSeq = —hVj is taken to be in the eigenstate 
representation, i.e. 7c = U^jU, with 
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(3) 



U as the k-dependent unitary matrix that diagonalizes 
the Hamiltonian at point k. 



(4) 



Sn is a matrix that is diagonal in the band indices, 
[Sn]ij = ^ij^im and the so-called Berry connection ma- 
trix is given by = iW d\JJ [11]. Not included in 
expression (|2]) are the vertex corrections, which vanish 



TABLE I: Anomalous Hall conductivities for bcc Fe and hep 
Co in units of S/cm for selected high-symmetry orientations 
of the magnetization, c^'''^ cr'"'^ and cr'^^^^''^ stand for IC, SJC 
and their sum, respectively. The experimental values are for 
the scattering-independent conductivity. 
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for an inversion-symmetric system in the Gaussian dis- 
order model. However, note that in contrast to the orig- 
inal formula as presented in Eq. (3) of Ref. [11], our 
expression for <j^j is manifestly antisymmetric. For the 
Rashba model it reduces to the original form of Eq. (3) 
in Ref. |ll]. It is important to note that the SJC in the 
short-range disorder model, Eq. (|2]), is solely determined 
by the electronic structure of the pristine crystal and thus 
directly accessible by ah initio methods. 

In practice, we replace the integrals in Eqs. ([T]) and ([2]) 
by a discrete sum over a finite number of /c-points in the 
Brillouin zone (BZ). To reduce the computational cost 
we adopt the method of Wannier interpolation [sl, [l^, 
which employs the description of the electronic structure 
in terms of maximally-localized Wannier functions (ML- 
WFs), to evaluate Eqs. ^ and (j2j) for bcc Fe, hep Co, 
fee Ni, as well as LIq FePd and FePt. The electronic 
structure calculations were performed with full-potential 
linearized augmented plane- wave method as implemented 
in the Jiilich DFT code FLEUR [HI within the generalized 
gradient approximation (GGA). We used the plane- wave 
cut-off iiTmax of 4.0 bohr~^ and 16000 /c-points for self- 
consistent calculations. Spin-orbit coupling was included 
in the calculations in second variation. We constructed a 
set of 18 MLWFs per atom using the Wannier 90 code [l6| 
and our interface between FLEUR and WannierQO [l7| . 

We present the results of our calculations of the in- 
trinsic and side-jump AHC for Fe, Co, Ni, FePd and 
FePt in Tables I, II, and Fig. [1] for high-symmetry di- 
rections of the magnetization M in the crystal. These 
results are compared to experimental values, from which 
the skew scattering contribution was either explicitly sub- 
stracted [qL llSf . or can be safely ignored at higher tem- 
peratures [i2[tli, [20| . 

We first analyze the results for bcc Fe. For M along 
[001], the IC in Fe accounts to roughly 75% of the known 
experimental value of ^^r^lOOO S/cm [6|, |20|. This compar- 
ison gets slightly improved considering that the experi- 
mental value averages over crystals with different orien- 
tation. Taking the SJC into consideration improves the 
value of the AHC in Fe significantly for all magnetization 
directions, with the angle- averaged cr^"* + of about 
90% of the experimental conductivity. In hep Co, the 
magnitude of the SJC for M along the c axis is as large 



as 217 S/cm, with the total AHC of 694 S/cm, very close 
to the experimental value of about 800 S/cm. For M in 
the basal ab plane the SJC is small and negative, bring- 
ing thus the intrinsic value down to ~70 S/cm, somewhat 
away from the experimental value of about 150 S/cm. 
One has to keep in mind, however, that the experimental 
values for hep Co are approximate [12 1. 

A significant improvement upon including the SJC is 
also evident for the more complex ordered FePd and FePt 
alloys in their LIq phase with M along the [001] axis. For 
FePd, the IC is very small, of about 130 S/cm, while the 
side-jump AHC is twice as large and of the same sign, 
resulting in a value of the total AHC much closer to ex- 
periment. Table II. As follows from our calculations, in 
FePd the AHC is dominated by a^^ , in accordance to an 
earlier indirect prediction [9|. On the other hand, in FePt 
with M|| [001], the IC is much larger, while the SJC is half 
the value of that in FePd. This is again in agreement to 
Ref. [qI, in which such a crossover between the intrin- 
sic and side-jump conductivities, appearing within the 
Dirac model as well [21|, was attributed to different SOI 
strength of Pd and Pt atoms. As far as the comparison 
to experiments is concerned, also in FePt adding the cal- 
culated SJC to the IC brings the total AHC within the 
range of experimentally observed values for samples of 
[001]-magnetized LIq FePt with high degree of ordering, 
S > 0.7, and different sample thickness P. fisj. 

The case of fee Ni presents a special challenge, since in 
this material the GGA value of the intrinsic AHC is much 
larger than the measured scattering-independent value, 
implying a sizable a^^ with the sign opposite to the IC ^ . 
From our calculations of the IC in fee Ni we obtain values 
which lie between -2000 and -2500 S/cm (Fig. [Hat U = 
0), depending on the direction of M, while the experi- 
mental value resides in the vicinity of —640 S/cm (il. [22j. 
The calculated values for the scattering-independent SJC 
in fee Ni presented in the same figure (at U = 0) lie be- 
tween — 100 and 400 S/cm, and thus cannot explain the 
large discrepancy between theory and experiment. The 
description of the electronic structure of Ni within con- 
ventional DFT is well-known to be inaccurate, however. 
Several attempts aiming at improving the GGA values 
for quantities such as magnetocrystalline anisotropy en- 
ergy, spin- wave dispersion etc., were made in the past 
(e.g. Refs. and references therein), proving the 
importance of correlation effects in this material and the 
sensitivity of calculated quantities on the shape of its 
Fermi surface. 

TABLE II: Same as in Table I for Llo FePd and FePt alloys. 
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FIG. 1: (Color online) The dependence of the intrinsic 
and side-jump (inset) conductivities on the strength of the 
Coulomb repulsion U of valence d-electrons in fee Ni for dif- 
ferent directions of the magnetization in the crystal. 



In our work, we choose the GGA-\-U approach in order 
to study the effect of correlations on the AHE in fee Ni, 
following the implementation of Ref. [25| and treating the 
double counting corrections within the atomic limit [26| . 
For this purpose, we scan the strength of intra-atomic re- 
pulsion parameter [/, keeping at the same time the intra- 
atomic exchange parameter J such that the value of the 
spin moment of Ni stays roughly constant ^| . As can be 
seen from our calculations, presented in Fig. 1, the values 
of the cr^"* upon including U change drastically and come 
closer to experiment, approaching a value of —800 S/cm 
when U is changed in the range of 0—4 eV, comm only 
used for calculations of other properties of Ni [IJ, \2q . 
This suggests that the main reason for the discrepancy 
between the intrinsic AHC values obtained from DFT 
and experiment might lie in the improper description of 
Ni's electronic structure from first principles. On the 
other hand, the values of the cr^^ are affected differently 
by the modifications in the electronic structure, almost 
not changing for M along the [001] axis, and displaying 
a non-monotonous behavior within the range of 100-300 
S/cm in the absolute value as a function of U for two 
other magnetization directions. 

Such different sensitivity to the band structure can be 
understood by analyzing the structure of the SJC and 
the IC on the Fermi surface (FS). To simplify things, 
when taking into account the very complicated FSs of 
the ferromagnets considered in this study, we sum up 
the contributions to a^"* and <j®j over all bands and all 
sheets of the FS, respectively, following Eqs. (1) and (2), 
when going along a certain direction in the BZ until the 
inner "Fermi" sphere in the BZ with the center at its 
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FIG. 2: (Color online) Angle-resolved conductivity da/dQ in 
units of S/cm as a function of direction in the BZ. da/dQ 
corresponds to all contributions to a from inside the inner 
sphere in the BZ within the solid angle element dQ. (a) a^^^ 



for Ni [001], (b) cr'J for Ni [001], 
for Fe [001]. 



(c) cr'J for Ni [110], (d) a'^ 



origin (F-point) is reached. In the case of such angle- 
resolved IC in Ni ([/ = 0), shown in Fig. 2(a), large 
contributions can be seen along the "hot loops" in the 
BZ, which are situated in the vicinity of the intersections 
between different sheets of bands [l3|, while the IC in 
the region away from such band crossings is also signifi- 
cant [sl . This is rather different from the topology of the 
SJC on the Fermi sphere. The SJC in Ni {U = 0) and Fe, 
presented in Fig. 2(b)-(d) manifests that the main con- 
tribution to cr^J comes from certain isolated " hot spots" , 
distributed rather sparsely over the Fermi sphere, while 
the SJC decays very quickly with the distance from such 
points. Such a strong difference in the distribution of 
the cr^j and a^"* on the FS arises from the effective mag- 
netic monopole nature of the IC Berry's phase contri- 
bution near the band crossings, resulting in a more pro- 
nounced sensitivity of the a^"* to the parameters of the 
electronic structure, such as Coulomb repulsion /7, Fermi 
energy etc., whereas the SJC does not contain such sin- 
gularities near those crossings. 

From our calculations presented above in Tables I, II, 
Figs. 1 and 2, it is evident that cr^^ exhibits large changes 
when the direction of the magnetization in the crystal 
is varied. In uniaxial crystals, such as FePt and hep Co 
such anisotropy is not surprising, given that in these ma- 
terials also the anisotropy of the intrinsic AHC appears 
already in the first order with respect to the directional 
cosines of the magnetization j^. And while the differ- 



ence in the absolute change in the SJC and IC in FePd 
and FePt upon rotating the magnetization direction can 
be probably related to different SOI strength of the two 
materials [9], in FePt the corresponding trend of the SJC 
and the IC is opposite owing to the different Fermi sur- 
face topology of the two contributions. Surprisingly, the 
strong anisotropy of the SJC can be also observed in bcc 
Fe and fee Ni. In Fe this anisotropy reaches as much 
as 70%, while in Ni the SJC anisotropy is striking as 
compared to the IC anisotropy, with changes in sign and 
order of magnitude of cr^^ as a function of the magneti- 
zation direction. This can be perhaps intuitively under- 
stood considering that a^^ is given almost entirely by sin- 
gular "hot spots" at the FS, which change their position 
and the magnitude of their contribution depending on 
the matrix elements of the SOI, controlled in turn by the 
magnetization direction [l2|, |28|, compare e.g. Fig. 2(b) 
and (c) for Ni. 

Despite the unprecedented improvement of the values 
of the AHC in several ferromagnets when compared to 
the experimentally measured numbers, the scattering- 
independent SJC considered here cannot describe the en- 
tire physics of the complex side-jump scattering and will 
likely fail to describe it in certain systems where long- 
range scattering and spin-dependent scattering domi- 
nate. This can be particularly important for the case 
of low-doped, i.e., having long screening length, magnetic 
systems. Also, within our approach, we consider only the 
leading-order in impurity strength correction to the self- 
energy, which is justified within the weak scattering limit. 
This approximation might fail, however, when the per- 
turbation in the crystal potential due to the presence of 
disorder or impurities is very strong. It would be highly 
desirable to extend the current model for the scattering- 
independent side-jump conductivity beyond the short- 
range disorder and weak scattering limit. 

In summary, we have implemented a method to cal- 
culate the scattering-independent SJC within DFT. We 
found that for fundamental ferromagnets, such as Fe, Co, 
FePd and FePt the agreement between theory and exper- 
iment can be essentially improved upon considering the 
scattering- independent SJC. This SJC can be calculated 
from the electronic structure of the pristine crystal only, 
which encourages the application of the considered model 
for the side-jump scattering to wider classes of materials 
with the goal of extending the applicability of the DFT 
in treating complex transverse scattering phenomena and 
comparison to experiments that are commonly performed 
on samples with unknown disorder and impurity content. 
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